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ABSTRACT 



We present Spitzer-IRS spectra obtained along the molecular jet from the 
Class source L1448-C (or L1448-mm). Atomic lines from the fundamental 
transitions of [Fell], [Sill] and [SI] have been detected showing, for the first time. 
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the presence of an embedded atomic jet at low excitation. Pure rotational H2 
lines are also detected, and a decrease of the atomic/molecular emission ratio 
is observed within 1' from the driving source. Additional ground based spectra 
(UKIRT/UIST) were obtained to further constrain the H2 excitation along the 
jet axis and, combined with the 0-0 lines, have been compared with bow-shock 
models. From the different line ratios, we find that the atomic gas is characterized 
by an electron density Ug ~ 200-1000 cm~^, a temperature Tg < 2500 K and an 
ionization fraction < 10^^; the excitation conditions of the atomic jet are thus 
very different from those found in more evolved Class I and Class II jets. We also 
infer that only a fraction (0.05-0.2) of Fe and Si is in gaseous form, indicating that 
dust still plays a major role in the depletion of refractory elements. A comparison 
with the SiO abundance recently derived in the jet from an analysis of several 
SiO sub-mm transitions, shows that the Si/SiO abundance ratio is ~100, and 
thus that most of the silicon released from grains by sputtering and grain-grain 
collisions remains in atomic form. Finally, estimates of the atomic and molecular 
mass flux rates have been derived: values of the order of ~10~^ and ~10~^ 
M0yr~^ are inferred from the [SI]25/im and H2 line luminosities, respectively. A 
comparison with the momentum flux of the CO molecular outflow suggests that 
the detected atomic jet has the power to drive the large scale outflow. 

Subject headings: stars: formation — ISM: jets and outflows — ISM: individual 
(L1448-C) — infrared: ISM — ISM: lines and bands 



Introduction 



The process of mass accretion, leading to the formation of solar type stars, is always 
associated with mass ejection in the form of collimated jets, which e xtend from a few AU up 



to parsecs away from t he exciting source. According to the models (IKonigl &: Pudritzl 12000 



Casse fc Ferreirall2000l ). accretion and ejection are intimately related through the presence 
of a magnetized accretion disk: the jets carry away the excess angular momentum, so that 
part of the disk material can move towards the star. This paradigm of star formation is 
now being observatio nally tested in Cl ass I and Class II objects through detailed optical and 
near-IR observations (IRay et al.l 120071 ). 



However, the characteristics of jets from these evolved YSOs are unlikely to be appropri- 
ate for those from protostars in earlier evolutionary phases, which are expected to propagate 
in a denser medium and be associated with more energetic mass ejection. In such un-evolved 
objects, so-called Class sources, the initial part of the jet is often detectable only at mm 
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wave l engths in the form of a collimated, hi gh- velocity molecular outflow (iGueth fc Guilloteau 
19991 : iLee et al.l 120071 : ICodella et al.l 120071 ) . Near-IR H2 emission is always obscured by high 
extinction near the central engine and jet base, and instead traces hot gas excited further 
downstream in bow shocks near or at the jet apex. While it is usually assumed that the 
molecular jet represents the cold external layer of an embedded atomic jet, in principle the 
jet could be composed entirely of low excitation molecular gas and this needs to be tested 
observationally. 

Millimeter interferometric observations of high velocity molecular jets from Class 
sources present characteristics very similar to the hot jets seen in T-Tauri or Class I YSOs, 
such as a very narrow (< 2'') width and a kn otty structure resembling that of HH objects 
( IGueth &: Guilloteaulll999l : ICabrit et al.l 120071 ). Sub-mm observat ions in CO and S iO show 
that such jets are dense, with peak values of n(H2) ~10^ cm~^ (e.g lNisini et al.ll2007l ). More- 
over, ISO-LWS observations suggest that they are also warm, with temper atures between 



300 and 1500 K inferred from the copious high-J CO and H2O emission (IGiannini et al. 



2OOII ). Such a warm gas component may represent the bulk of the mass flux ejected by the 
protostar, and thus, energetically, may be the jet's most important component. The very 
low LWS spatial resolution, however, has not permitted us to draw any conclusions about 
the detailed structure of this warm gas, giving only physical parameters averaged over the 
entire outflow. 

Spitzer observations now allow us to investigate the properties of the warm gas com- 
ponent of the molecular jet through mid-IR atomic and molecular features probing low 
excitation gas at T~100-2000 K and n ~10^-10^ cm~^. The moderate IRS spatial resolution 
(4-20") allows one to separate, in nearby sources, the inner jet region from the region where 
the jet is strongly interacting with the ambient medium through bow-shocks. 

In this study, we present results of an analysis performed on Spitzer-IRS spectra of 
the jet from the Clas s source L1448-C (or L1448-mm). This is a low luminosity (L=7.5 
Lq, iTobin et al.l 120071) pro tostellar source located in the Perseus Molecular Cloud (D~ 250 
pc, e.g. Enoch et al.ll2006l ). A powerful and highly c ollimated flow is driven by this source. 



as te stifled by interferometric CO and SiO maps (iGuilloteau et al.l Il992l : iBachiller et al. 



I995I ). This highly collimated molecular jet is associated with a less collimated and ener- 
getic CO outflow, probably representing ambient swept-up material. Near-IR observations 
performed on the L1448-C outflow show that H2 hot molecular gas is detected only at the 
bow shocks at a distance of ~1' from the central source; the underlying jet that is responsible 
for the shocks remains undetected at optical and near-IR wavelengths a nd has been traced 
at mm- wavelengths as collimated SiO emission (IGuilloteau et al.lll992l ). Far-IR observa- 
tions performed with ISO have shown the existence of a warm gas component (T~1000 K, 
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n~10^ cm ^) associated with the molecular jet, s howing; up in [OI]63/im, H2, CO and H2O 



pure rotational emission (INisini et al.lll999l . l2000l ). Recent Spitzer IRAC ar id MIPS images 
have shown that t he so urce is actually binary, with a separation of ~7" (|j0rgensen et al. 

20071). The north/south binary components are called L1448-CN and 



20061: 



Tobin et al.l 



CS (|j0rgensen et al. 2006) or L 1448-mm A and B ( Tobin et al. 2007) respectively. We will 



use the |j0rgensen et al.l (120061 ) nomenclature in this paper. Of the two sources, the CN 
component is associated with the strong millimeter source rec ognized by interferometric 
observations as the outflow driving source (IBachiller et al.lll995l ). 



The paper is organized as follows: §2 describes the Spitzer observations and data reduc- 
tion as well as additional near-IR observations obtained at the UKIRT telescope. §3 and 4 
contain our analysis of the observed H2 and fine structure lines in order to derive the main 
physical parameters of the jet. Finally, conclusions are presented in §5 . 



Observations and data reduction 



2.1. SPITZER - IRS 



The driving source and outflow of the Class object L1448-C were observed with the 
SPITZER Infrared Spectrograph (Houck et al. 2004) during March and September 2006. 
Four positions have been targeted with the Short-High (SH), and Long-High (LH) IRS 
modules [R ~ 600, 10-37yum) in order to cover the central driving source (L1448-CN), the 
southern CS source (the position of which spatially coincide s with a CO red-shif ted clump 
of Extremely High Velocity (EHV) jet emission named Rl (IBachiller et al.lll990l ). and two 
adjacent positions along the blu e-shifted part of the je t, that we call OFl and 0F2, the first 
comprising the Bl mm clump in Bachiller et al. Il990) and the sec ond the H2 bow-shock, as 



detected in the near-IR image of the flow by lDavis fc Smith! (119961 ). respectively (see Figure 
d]). The plate scale of the detector is 2'.'3 pixel"^ and 4'.'5 pixel"^ for the SH and LH modules, 
respectively. Observations with the Short-Low (SL) module {R ~ 60-120, pixel scale 1'.'8 
pixel"^) were also acquired to extend the spectral coverage down to 5/im. The observing 
period was chosen in order to have the SL slit aligned along the L1448 jet. The observations 
were performed in staring mode with a total integration time of 2.5 hrs per position. 

The data were reduced and calibrated with the S13 pipeline. Basic calibrated data (bed 
images) were median-combined and cleaned for rogue pixels with the IRSCLEAN_MASK 
routine, and additional bad pixels in the low resolution module were removed by visual in- 
spection. Spectral extraction was performed with the Spitzer IRS Custom Extraction tool 
(SPICE). For the high resolution modules the full slit was extracted, whereas for the SL mod- 
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ule the full slit length was divided into four equal regions (~ 14") that were consecutively 
extracted in order to best match the high resolution module pointings. High resolution spec- 
tra were defringed using the IRSFRINGE package while the zodiacal light contribution was 
estimated using the Spitzer Planning Observations Tool (SPOT) and subtracted. Inter-order 
flux offset and curvature effects were minimized by optimally selecting the best approxima- 
tion between the point and extended source calibration options (SPICE), according to the 
morphology of each observed region. The resulting line flux difference between the two ex- 
traction methods is ~ 20% which can be considered as an upper limit to the error for the 
observed line fluxes. 

The calibrated spectra combined for all modules at each position are presented in Figure 
[21 in the same figure open squares and triangles represent IRAC bands 3,4 and MIPS band 



1 cont inuu m flux measureme nts from the recent studies of the region by iJcirgensen et al. 



(120061 ) and lTobin et al.l (120071 ) respectively. Their good agreement with the extracted spectra 
baseline confirms the appropriateness of the adopted extraction and calibration techniques. 
In particular, the procedure for zodiacal light subtraction is sufficiently precise to cancel the 
baseline jumps between different modules, the only exception being the SH and LH modules 
at the CS position, around 20/im. 

The extracted spectra display a number of emission lines that arise both from molecular 
and forbidden atomic or ionic transitions. Line identification has been performed considering 
the features lying more than 3a above the local rms and whose emission peak, fitted by a 
gaussian, is within half a resolution element from the theoretical vacuum wavelength. The 
observed molecular emission lines are due to H2, where the full series of pure rotational 
transitions (S(0) - S(7)) in the given wavelength range of the instrument are detected. Atomic 
and ionic emission has been detected in the form of the fundamental fine-structure lines from 
[Fell] (25.98/im), [SI](25.25/im) and [Sill] 34.81/im (see Fig|3D. None of the other numerous 
higher excitation [Fell] lines falling in the investigated wavelength range has been detected 
at the 3a level: this fact, combined with the non-detection of the [Nell] line at 12.8 fim, 
already suggests that the investigated region presents very low excitation conditions. 

In addition to emission lines, wide absorption bands are also observed in the spectra 
extracted at the two sources, the most prominent being those due to silicates at 9.7 fim, water 
ice at 6.02/im and CO2 ice at 15.2/im. 

Fluxes of the individual lines were computed by gaussian fitting after subtracting a 
sloped baseline with in the Spectroscopic Modeling, Analysis and Reduction Tool (SMART) 



(IHigdon et al.ll2004l ). In Tabled] we report the parameters of the detected lines, along with 
their identified transitions and the associated rms errors given by the la uncertainty derived 
from the fluctuations of the local baseline. The measured fluxes have a larger uncertainty of 
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up to 20% given by the flux calibration uncertainty. The 0-0 S(2) line has been observed by 
both the SL and SH modules. The line fluxes measured in the two modules agree inside the 
errors with the exception of the position 0F2, where the SH module gives a flux higher by a 
factor of two. This difference could be due to a different bow-shock area sampled by the two 
modules. For consistency, we use in all the positions the average of the two measurements. 

In the same table we also give the upper limits for a set of non-detected lines that fall 
within the observed range and that will b e used in the analysis that follows. In a previous 



study of the same region with ISO/SWS, iNisini et al.l (119991 ) detected the S(3) to S(5) H2 



transitions in a region of 14" x 27" centered on the L1448-C source. The derived fluxes 
are a factor between 3 and 8 larger than the ones presented here. The FOV of the SWS 
observations is around 10 times larger than the region extracted in our spectra: therefore 
the observed flux difference can be attributed to the different adopted apertures. 



2.2. UKIRT - UIST 

In order to further investigate the excitation conditions along the jet axis, additional 
spectra, covering the wavelength range from 1.4 to 2.5/im, were obtained on 4th October 
2006 at UKIRT (UK Infrared Telescope) using the image-spectrometer UIST in spectro- 
scopic mode. A 4-pixel-wide slit was used with a pixel scale of 0.12", corresponding to a 
spectral resolution ~ 500-700. The slit was positioned roughly parallel with the jet axis and 
thus aligned with the SL IRS slit, at a position angle (PA) of -17°. Object-sky-sky-object 
sequences were made with a total exposure time of 300s for each one. Each spectral image 
was bias subtracted and flat-flelded using the ORAC data reduction software. Further re- 
duction was performed using IRAF standard tasks. An Argon lamp was used in order to 
wavelength calibrate the spectra. B-type stars were observed with the same configuration in 
order to remove the telluric features and flux calibrate the spectra. Once the reduction was 
completed, the APALL IRAF task was used to extract the spectra corresponding in length 
to the four regions investigated with IRS. Table [2] lists the lines detected in the CN, OFl 
and 0F2 regions, with the corresponding measured flux. No lines have been detected in 
the CS position. Only H2 1-0 and 2-1 ro- vibrational lines were detected, while any atomic 
or ionic emis sion (e.g. the strong [Fell] lin es at 1.25/xm and 1.64/im) is missing, as already 



evidenced by lCaratti o Garatti et al.l (120061 ). It can be noted that in the 0F2 position, which 
corresponds to the bow-shock position, a large number of transitions from higher excitation 
1-0 lines (up to S(9)) are detected. This may be indicative of higher excitation conditions 
associated with the bow-shock or lower reddening pertaining to this region since it is located 
far from the mm source core. 
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3. Analysis of the excitation conditions along the jet 
3.1. Extinction and spatial variation of line luminosity 



Optical extinction (Ay) values along the line of sight of the CS and CN sources have 
been measured from the 9.7/im silicate absorption feature, adopting the relationship be - 
tween the optical depth of this feature and the visual extinction given by iMathid (119981 ) . 
i.e. Av/Tg/i=19.3 mag. The derived Ay values are 32 and 11 mag for the CS and CN 
source, respectively. In the OFl and 0F2 positions no silicate absorption is detected so 
no direct measurement of Ay fro m our spectra is po ssible. We have therefore assumed an 
extinction of 5 mag, as derived by lNisini et al.l (120001 ) towards the direction of the Bl clump 
(coincident with the OFl position). Such a low value of extinction is in agreement with the 
absence of water and CO2 ice absorption features in the spectra of these positions, as these 
features should become detectable in the Spitzer spect ra only once the visual extinction 
reaches values above 4.3±1.0 mag (IWhittet et al.l 120071 ). All the individual detected lines 
have been dereddened adopting these Ay values and the different extinction laws appropriate 
for the considered wavel ength ranges: in the range from 1 to 13 fim the extinctio n law of 
Rieke fc LebofskyI (119851 ) was adopted while from 13 to 23 fim the extinction law of iMathis 
( 119981 ) was used, which was extrapolated to 28 fim. 



An aposteriori check of the correctness of the adopted extinction values was done using 
the Boltzmann diagrams constructed from the H2 near and mid-IR lines (see Sect. 3.2), by 
examining the fit to a straight line of transitions appropriate for the 9.7/im silicate absorption 
(like the 0-0 S(3) line at 9.67/im) and at the interface between the near- and mid-lR lines 
in the same plot. 

Once extinction corrected, we have explored the spatial variation of the relative atomic/molecular 
emission, by plotting the ratio of the [Fell]26/im and [Sl]25/im with respect to the H2 0-0 S(l) 
line (Figured]). This plot clearly shows that the relative brightness of the atomic component 
with respect to the molecular component sharply decreases going from the CN to the OFl 
position while only slightly increasing again at the bow shock position. A similar behavior 
has been observed in the near-lR for a number of Class 1 jets, where the atomic and molec- 
ular components have b een traced by the [Fell] 1.64/xm and H2 1-0 S(l) lines respectively 
( iNisini et al.ll2002l . |2005| ). In these sources, the relative decrease of atomic gas emission with 
respect to H2 in the jet beams is accompanied by a decrease of excitation in the jet which 
occurs on scales of ~ 5-10" from the driving source. The interpretation of this behavior is 
that the jet, heated and ionized in the acceleration region, progressively expands and cools 
down until it strongly interacts with the medium through a bow shock. The spatial scale of 
the IRS instrument is too poor to resolve the intensity and excitation variations within the 
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jet expansion region; a qualitative inspection, however, suggests a similar behavior, in spite 
of the different physical conditions pertaining to the L1448 jet. 

3.2. H2 emission 

Being very easily thermalized, H2 lines can be used as probes for the temperature by 
means of excitation diagrams. The latter are constructed by plotting the values of ln(N^^j/g) 
against E^^j, where is the column density for the population at the upper level, E^^j and 
g are the level energy and statistical weight. Assuming LTE conditions, these two quantities 
are linearly related, and the local temperature can be derived from the slope of the fit. We 
have constructed excitation diagrams in each position combining the Spitzer and UKIRT 
observed lines. For these, column densities have been derived dividing the line flux for the 
aperture adopted for the spectral extraction, assuming that the emission is extended and 
fills the extraction apertures in a uniform way. This cannot be the case for the 0-0 S(0) line 
observed with the much larger aperture of the LH module: the H2 fundamental transition 
should be in fact also affected by diffuse emission from the cloud, therefore the derived 
column density can be considered as an upper limit of the S(0) column density in the jet. 
In constructing the Boltzmann diagrams, we have assumed an H2 ortho/para ratio equal to 
the equilibrium value of 3. We do not appreciate significant variations from this value inside 
the errors of the data-points. 

The excitation diagrams for each position are presented in Figure [5l in these, Spitzer 
(open circles) and UKIRT lines (open squares) are least square fitted with a straight line. 
The visual extinction values used and temperatures derived from the slope of the fitted lines 
are displayed in the top right of each diagram. In all cases, data points are reasonably 
well aligned and thus giving support to the assumption of a H2 gas in LTE. Because of the 
high values of visual extinction at CS, the lack of data points results in a fit with large 
uncertainties in the derived temperature. 

Given the fact that the lines falling within the Spitzer and UKIRT ranges originate from 
pure rotational and ro-vibrational transitions, respectively, the two sets of lines may trace 
different temperature regimes. However it is only in the bowshock position that the two 
sets of lines clearly trace two different temperature components. In the other positions, the 
extinction corrected upper limits on high excitation H2 lines suggest that a single temperature 
component with T ranging between 600 and 900 K, well represents the excitation conditions. 
It is only in the OFl position that the column density of the S(6) and S(7) lines deviate 
from the single temperature fit: we think that the displacement of these lines is not reliable, 
since such curvature is not followed by the NIR measurements of the higher excitation lines. 
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From the Boltzmann diagrams it is also possible to derive the total H2 column density of the 
warm gas, from the intersection of the fitted line to the data points and the abscissa of the 
excitation diagram, which is equal to InN {H2) / Q{T) , with Q{T) is the partition function at 
temperature T. 

The derived temperatures and column densities for each of the observed regions are 
presented in Table [31 The gas temperature shows a slight decrease going from the driving 
source (CN) to the OFl position, while column density displays a constant increase from 
the CS to the 0F2 positions. This trend confirms that the excitation conditions decrease as 
the jet propagates into the ambient medium, and that in the inner region the jet is mostly 
atomic, as testified from the [FeII]/H2 ratio displayed in Figure HI 



3.2.1. Shock model 

As seen in the previous section, the H2 Boltzman diagram in the 0F2 position shows 
the presence of different temperature components, likely arising in the unresolved cooling 
zone behind the bow shock. In order to constrain the shock conditions from our Spitzer and 
UKIRT observations, w e have attem pted to model the derived H2 column densities with a C- 



type bow-shock model (ISmithlll99ll ). For a better comparison of the derived column density 
with the model fit, we have constructed a Column Density Ratio (CDR) plot, in which the 
column density in the upper energy level of each transition has been normalized to the value 
given by a gas at lOOOK. These values are then presented relative to the column of the 1-0 
S(l) upper energy of 6953K (see Fig. [6]). In this Figure, triangles represent Spitzer 0-0 lines, 
squares represent all lines originating from the first vibrational levels, including the 1-1 upper 
limits derived from the Spitzer IRS spectra, and diamonds correspond to 2-1 lines; model 
predictions are superimposed over the data. The predicted column densities are calculated 
assuming that the upper energy levels are in rotational LTE but vibrationally in non-LTE. 
This is a reasonable assumption for the levels at low excitation energy, while the higher 
rotational levels appear to switch across from the 1-0 line to the 2-1 values, indicating that 
this assumption probably does not hold anymore for the higher levels. The best fit model 
prediction is obtained with a bow shock moving at 100 km s~^ into a medium of (H nucleon) 
density 10^ cm~^. The ion fraction of 10~^ and magnetic field strength, corresponding to an 
Alfven speed of 2 km s~^, imply that the non-dissociative wings of the bow are C-type shocks. 
The bow geometry required for the obtained fit is such that the bow must possess a compact 
apex and extended flanks with a shape Z (x W with s = 1.38 in cylindrical coordinates. 
Then, a large column of low excitation gas is generated by the bow. A C-bow shock model 
with similar parameters as adopted here was able to reproduce also the kinematic and the 
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1-0 S(l) morphology of the considered bow shock (jPavis fc Smith! Il996l ). 



3.3. Atomic Jet 

As pointed out in Section 12.11 the detection of only the fundamental [Fell] transition 
indicates that the excitation conditions of the atomic jet are very low. The [Fell] 26 /im line 
originates from the level, that has an excitation temperature of ~550 K, while the 

level just above this one, gives origin to a transition at 35 fim, having T^x ~ 960 K. In 

principle, the non detection of the 35 /im transition could give a strong constraint on the gas 
temperature. However, the spectral region around 30 yum is very noisy and the upper limits 
are not stringent. To get constraints on the temperature, we have instead considered the 
upper limit on the ^^7/2 — ^ transition at 18 /im, that lies in a region of lower noise. In 
addition, we have also used the ratio [SiII]35 /im/[FeII]26 /im as a density probe. Si and Fe 
have a comparable ionization potential (8.15 and 7.9 eV respectively) thus a similar degree 
of ionization. Furthermore these two lines are excited at a similar temperature of ~ 500 K, 
thus their ratio depends only on the electron density and on the [Si/F e] gas phase abundanc e 



ratio. We assume a [Si/Fe]c,as ratio equal to the solar ratio (taken from lGirart fc Acordll200ll ). 
an assumption implying that the two species are equally depleted on grains. The validity of 
this assumption is further examined in Section HI 

Figure [7] presents a plot of the [SiII]34.8//m/[FeII]26.0/im ratio as a function of the 
electron density, for temperatures of lOOOK and 2000K (solid and dashed lines). This diagram 
has been constr ucted employing a statistical equilibrium model that considers the first 16 



levels for [Fell] ( iNisini et al.ll2002l ). and a two level system fo r [Sillj . Radiative and collisional 



rates for the [Sill] have been taken from lDufton fc Kingston! (jl99ll ). Hatched areas represent 



the ratios observed in the CN and 0F2 regions, respectively. Such observations are consistent 
with values of Ue in the range ^^300-500 for the 0F2 position and ~200-1000 cm~^ for the 
CN position. 

In Figure El the upper limit ratio [Fell] 18/im to 26/im is plotted as a function of the 
gas temperature, for different electron density values. The 26/im line flux and the 18/im 
line upper limit have been measured on two IRS modules having different field of views 
and therefore their intrinsic ratio depends upon the extension of the emitting region. A 
conservative upper limit on this ratio is obtained assuming beam filling and consequently 
normalizing the line ratio to the different field of view. The upper limits obtained in this 
way are displayed in Fig. [8] for the CN and 0F2 regions. From this plot, we derive that gas 
temperatures less than ~ 2500 K and ~ 1500 K are responsible for the emission observed in 
the 0F2 and CN positions, respectively. 
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The derived physical conditions significantly differ from the conditions measured from 
[Fell] near-IR lines in jets from more evolved class I sources, that have temperatures ranging 



from 7000 and 15000 K and densities - 10^ -10^ cm'^ (e.g. iNisini et all 12002 : IXakami et al. 



20041 ) . Low electron density values, in particular, point to a low ionization fraction for the 
gas under consideration. Total nH2 densities of the order of 10^ cm~^ or higher have bee n 
inferred in the L1448 jet from sub- mm and far-IR observations (INisini et al.l l2000l . 120071 ). 
implying a ionization fraction Xg of ~ 10~^ or lower. 

The inferred low temperature and ionization fraction, in conjunction with the non- 
detection of the [Nell] 12.8 /xm line, imply that, if these lines are shock excited, the shock 
velocity should be low. Fast J-type shocks (e.g. HoUenbach & Mc Kee 1989) are therefore 
excluded. Low-velocity J-type shocks might provide sufficient ionization to excite the funda- 
mental ionic fine-structure lines, but the shock velocity needs to be lower than ~ 10 kms~^ 
to have a temperature of less than ~2000 K. 

Part or all of the atomic emission observed on-source may in principle also originate 
from excitation in a circumstellar disk. Line emission disk models (Gorti & HoUenbach 
2008) as well as Spitzer observations in T Tauri disks (e.g. Pascucci et al. 2007, Lahuis et 
al. 2007) show however that disk emission is characterized by both strong [Nell] 12.8 fim, 
excited in the high temperature gas heated and ionized at the disk surface, and [Fel] 24 fim 
line, originating in the deeper disk vertical layers, that we do not detect here. 

The strong atomic emission observed on-source may instead originate from the jet base, 
in a zone similar to the Forbidden Emission Line (FEL ) regions observed in Class I/II jets 
although with very different excitation conditions (e.g. iDavis et al.ll2003l ). In analogy with 
the FEL regions there is a decrease in intensity of the atomic emission with distance from 
the source, likely caused by the expansion of the jet and the consequent cooling down of the 
gas. 



4. Dust disruption and abundance of refractory species 



Gas phase abundances of Iron and Silicon are very low in the interstellar medium since 
these refractory elements are easily depleted onto the cores of dust grains. Shocks occurring 
along the outflows of young stars are able to at least partially restore the refrac tory element s 
to the gas phase, through processes like sputtering and grain-grain collisions J JonesI booof) 



Observations of mm and sub-mm SiO lin es in many molecular outflows (e.g. ICodella et al. 



20071 : iGibb et al.l 120071 : INisini et al.l 120071 ) have shown that dust grains are indeed partially 



destroyed, and that the released Si undergoes chemical reactions leading to the formation 
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of SiO. flSchilke et al.l 119971 : iGusdorf et al.ll2008l ). SiO abundance determinations along the 
molecular jet of L1448 show however that gas-phase Si locked in SiO is only about 5 10~^ 
of the Si solar abundance: thus either not all of the Si released by the dust reprocessing is 
converted into SiO, or the shocks are not able to completely restore all the Si to the gas 
phase. The detection of the [Sill] fundamental line in our spectra suggests that indeed a 
significant part of gas-phase Si is present in ionic form. On the other hand, the gas-phase 
abundance of Fe in jets have been so far measured only in near-IR jets of C lass I sources, 
from the bright near-IR [Fell] lines ( iNisini et al.ll2002l . l2005l : iPodio et al.ll2006l ). Such studies 
have shown that a large fraction of Fe (from 70 to 95%) is still locked in grains, indicating 
that dust grains have not been totally destroyed by shocks. No estimates have been so 
far given about the iron gas-phase abundance in Class molecular outflows. We can now 
provide abundance estimates of both Si and Fe through the detected emission lines of their 
single-ionized atoms. 

The gas phase abundance of refractory elements can be derived from a comparison of 
their emission lines with those of a non-refractory species emitted under the same physical 
conditions. In the case of our spectra, the [SI]25.2/im line can be used as a reference, as 
sulfur is not depleted in grains, assuming that all the sulfur is in neutral form. To check 
this hypothesis, we have computed the S^/S"*" ratio applying ionization equilibrium between 
coUisional and charge-exchange ionization, and direct/ dielectronic recombination (rates from 
Standi et al.lll998l : iLandini fc Monsignori Fossil 1 19901 ). For the inferred physical conditions 
of T ~ 1000-2000 K and Xe < 10~^, it results that 95% of S is in neutral form. The iron gas 
phase abundance, relative to the solar value, can be therefore written as: 



[Fe/S] F([FeJJ]26/im) e([5/]25.2/im) 







[Fe/S]Q F([5/]25.2/im) e([Fe//]26/im) 



X [S/Fe] 







where e([S'/]25.2/xm) and e(\F eII]2Qum) are t he th eoretical emissivities, [Fe/ S]q is the 
solar abundance ratio, taken from lAsplund et al.l (l2005l ). and F([FeII]26/im)/F([SI]25.2/im) 



is the observed ratio. A similar expression can be written for the [Si\gas/[Si\Q ratio. To 
compute the theoretical emissivities, we have considered for S a 5 level statistical equilibrium 
code assuming electronic collisional excitation and ne=400 cm~^. Two values of temperatures 
have been considered: 2500 K (the upper limit derived in the 0F2 position) and 600 K (the 
lower value derived by the H2 analysis). In a medium with high total density and low 
ionization fraction, such as the one we are considering here, collisions with atomic hydrogen 
may become important in the excitation of atomic species such as [SI] , in spite of the low nn 
collisional de-excitatio n rates with respect to t he ra tes for electronic collisions {^uhe ~ 10~^ 
for the [SI]25/im line, iHoUenbach fc McKed ( 119891 )). In order to asses how collisions with 
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hydrogen affect the resuhs, we have calculated the emissivities also assuming a medium with 
a total density of 10^ cm"^ and uh = 0.1 Utot • 

We have applied the above analysis to the values observed in the CN, OFl and 0F2 
positions: for the [SI]25.2/im flux, we have taken the average value given by the determi- 
nations obtained in the two different orders of the LH module. Results are summarized in 
Table HI the uncertainty in temperature results in a factor of two the uncertainty in derived 
abundances, while a difference of four is found between results obtained assuming collisions 
with electrons and with atomic hydrogen. 

The Table shows that the gas phase abundance of Fe and Si remains between 5 and 20% 
of the solar values for both species and in all positions. Similar low values of Fe abundances 
have been also foun d in the inner regions of jets from C lass I sources, from the analysis of 
the [Fell] NIR lines (INisini et al.ll2005l : iPodio et al.ll2006l ). while larger abundances (from 20 
to 70% of the s olar value) have been derived at large d istances from the driving source and 
in bow shocks (INisini et al.ll2002l : iGiannini et al.l 120081 ). 



Finally, taking a SiO abundance of 5 10~^ derived by iGusdorf et al.l (120081 ) . the Si/SiO 
abundance ratio is ~ 100, giving support to the hypothesis that most of the silicon released 
from grains remains in atomic form. 



5. Mass and mass flux in the warm gas 



The presented observations point to the presence of warm gas associated with the L1448 
CO and SiO millimeter jet, composed of both a molecular and a weakly ionized component. 
We can assess if this warm gas represents a dynamically important component of the jet, 
by measuring its mass flux and comparing it with estimates of the mass flux based on ISO 
and sub-mm observations of the CO emission (INisini et al.ll2000l : iBachiller et al.lll990l ). We 
can determine the mass flux for both the molecular and atomic counterparts of the outflow, 
using as tracers H2 and [Fell] /[SI] respectively. The method is similar in both cases and is 
based on the fact that the line emission is opt ically thin, so that th e observed luminosity is 
proportional to the mass of the emitt ing gas (|Hartigan et al.lll994j ). For [Fell] and [SI] we 
have applied the relationship given in INisini et al.l (120051 ) : 



^neutral hydrogen collisional rates for ions, such as [Fell] and [Sill], are about 10^ times weaker than 
electron colhsional rates: we have checked that nn colhsions start to give comparable ion emissivities with 
respect to electron collisions only for nn > 10^ cm~'^ and Ug < 500 cm"'^. We have therefore not considered 
collisions with n^ in our calculations of [Fell] and [Sill] emissivities 
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M = finiH X (uhV) X (dvt/dlt) (2) 

where /i is the mean atomic weight, mii^ is the proton mass, njy is the total number 
density, V the volume of the emitting region, dlt is the projected length perpendicular to the 
line of sight and dvt the tangential velocity of the observed region. The number of emitting 
atoms can be derived from the line luminosity according to the relation: 



uh V = L{line) I huAif, 



X 
H 



(3) 



where Ai and ft are the radiative rate and fractional population of the upper level, and 
[X/H] is the gas phase abundance of the considered atom/ion. For the determination of the 
mass flux from the [FeII]26/im line, we have assumed that all iron is singly ionized based 
on the fact that no neutral iron lines are observed within the IRS range; this assumption is 
based on the fact that the [Fe I] ground transition at 24/im is not detected in any of our 
spectra, in spite of its high radiative rate coefficient, comparable to that of the [Fell] ground 
transition. We have taken the iron gas phase abundance estimated in Sect 3.1 and listed in 
Table HI Conversely, to apply the relationship (3) to the [SI]25/im line, we have assumed 
that sulfur is all neutral (as discussed in Sect. 3.1) and taken the sulfur solar abundance. 

The projected length dlt was taken equal to the jet length sampled by the width of 
the LH slit. Finally, the tangential velocity of th e outflow was taken eq ual to 170 km s~^ 



as derived from the SiO proper motion study by iGirart fc AcordI (120011 ). Such a velocity, 
appropriate for the sub-mm jet, may be just a lower limit for the atomic jet velocity, if 
the latter is associated with a higher velocity inner component. It has been shown however 
that atomic jets from Class I/II sources have velocities in the range 100-300 kms~^ (Ray et 
al. 2007, Davis et al. 2003). We therefore believe that our assumption can introduce an 
uncertainty of at most a factor of two in the mass flux determination. The derived mass 
flux values using the above method are listed in Table [H The uncertainty associated with 
the Fe gas-phase abundance results in an order of magnitude uncertainty in the mass flux 
determination using the [FeII]26/im. The determinations from the [SI] line agree with the 
range of values derived from [Fell] and indicate a mass flux of the order of 1-2 10~^ Mq yr~^ 
in the CN position, while lower values, of the order of 3-8 10~^ Mq yr~^, are estimated at 
the bow-shock 0F2 position. 



The determination of the mass flux from the H2 component has been measured in- 
stead from the total column density for each extracted region as derived from the excitation 



- 15 - 



diagrams, applying the relationship: 

M = irniH X {2N{H2)A) x (dvt/dlt) 



(4) 



where N(H2) is the total column density and A is the area sampled by the slit. We con- 
sider here again the tangential velocity as measured from the SiO proper motion, assuming 
that it is representative of all the molecular gas. The derived mass flux along the outflow 
with this method is presented in Table O The H2 mass flux in the central position is about 
two orders of magnitude lower than the value derived from the atomic lines while in the 0F2 
position M(H2) is a factor between 2 and 5 lower than M([SI]).Thus the bulk of the flowing 
mass is carried out by the atomic jet all along the sampled positions. 

// 

A mass flux of the order of 5 10^ MqJt~ has been der ived within 20 o f the driving 
source combining multi-transition CO sub-mm and ISO data (INisini et al.ll2000l ). This value 
is only a factor of few larger than the value we derive from atomic lines at the CN position, 
which suggests that the atomic jet detected by our observations and the warm CO outflows 
sampled by the ISO and sub-mm observations are dynamically linked. It is also instructive to 
examine if the atomic jet detected by our Spitzer observations posse ses enough momenturn 
to sustain the entrained CO outflow observed at mm wavelength. iBachiller et al.l (119901 ) 
derived, for the L1448-mm blue-shifted lobe, a momentum flux of the order of 2.6 10~^ 
Mq kms^^ yr^^, assuming an outflow inclination of 70° with respect to the plane of the sky. 
If we assume momentum conservation for entrainment of the outflow by the jet, and a total 
jet velocity of 170 kms~^, we derive that the primary jet should possess a mass flux of the 
order of 1.4 10~® Mq yr~^ . This is consistent with the values we infer from the atomic 
emission, supporting the hypothesis that the detected atomic jet represents the energetically 
most important component of the outflow. 

The accretion rate can be also estimated assuming the bolometric lumino sity is mainly 
due to accretion. Using a stellar mass ~ 0.5 Mq (IFroebrich et al.l 120031 ) and a stellar 
radius i?^ ~ 4 Rq (IStahler et al.lll980l ). we get an accretion rate Mace ~ Lhoi R-^/G ~ 
2 10~^ M0yr~^. We thus derive a value comparable to the mass ejection rate estimated 
above. This is inconsistent with the picture that about only 10% of the mass accreting 
on the protostar is ejected by means of magneto-centrifugal acceleration mechanisms (e.g. 
Ferreira et al.l 120061 ). This may indicate that the actual protostellar mass is lower than 
assumed, and that the source is still at the beginning of its main accretion phase. 
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6. Conclusions 

We have carried out Spitzer IRS spectroscopic observations towards the molecular jet 
driven by the Class source L1448-C. We have detected the H2 pure rotational lines from 
S(0) to S(7) alongside fundamental fine-structure transitions of Si+, S° and Fe+. Additional 
UKIRT/UIST spectra have been acquired, showing that, at variance with the 0-0 lines, 
near-IR H2 vibrational line emission is mainly confined to the bow shock region where the 
jet and the ambient molecular cloud interact. We have analysed these lines in order to derive 
the physical and dynamical parameters of the associated warm gas. The main results can 
be summarised as follow: 

• The detection of the fine-structure lines testify for the presence to a previously un- 
known atomic jet, embedded in the molecular outflow. An analysis based on the line 
ratios as well as on previous observations at sub-mm wavelengths indicates that the 
excitation conditions of this jet are very low. In the region close to the driving source, 
a temperature Te < 1500 K and an electron density ~200-1000 cm~^ have been 
measured, while the estimated ionization fraction is surprisingly < 10~^. We suggest 
that the detected atomic gas represents the analog of the Forbidden Emission Line 
(PEL) region observed in the more evolved Class I/II sources, although with very dif- 
ferent excitation conditions. The close similarity with the PEL region is seen in the 
decrease of intensity of the atomic emission with distance from the source; in both 
cases this is likely caused by the expansion of the jet and the consequent gas cool 
down. Observations with better spatial and spectral resolution, hke those that will 
be possible with the Mid-Infrared Instrument (MIRI) instrument on the James Webb 
Space Telescope (JWST), are needed to have a better picture of the origin of this gas 
through comparison with jet acceleration models. 

• The H2 rotational emission indicates the presence of warm gas at a temperature ranging 
between 600 and 900 K. It is only at the bow shock position that a second component 
at higher temperature is detected in the near-IR. The H2 emission in the inner jet 
region may represent warm molecular gas enveloping the atomic gas: such a molecular 
component can be also related with the coUimated SiO jet observed at sub-mm wave- 
lengths, the excitation conditions of which are similar to those inferred here for the 
H2. 

At the bow shock position, the near and mid-IR H2 line observations have been com- 
pared with a C-type bow shock model. Shock velocities of ~ 100 kms~^ and pre-shock 
densities n(H2) ~ 10^ cm~^ can account rather well for the observed column density. 

• The gas phase abundance of Fe and Si has been estimated to infer the amount of 
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dust reprocessing in the jet. Only a fraction of ~0. 1-0.3 of the solar abundance of Fe 
and Si is in gaseous form indicating that dust still plays a major role in the depletion 
of refractory elements. A comparison with the SiO abundance recently derived from 
the analysis of sub-mm transitions shows that the Si/SiO abundance ratio is ~100, 
and thus that most of the silicon released from grains by sputtering and grain-grain 
collisions remains in atomic form. The inferred Fe gas-phase abundance is similar to 
the values previously estimated in Class I jets by means of near-IR lines: this suggests 
that dust reprocessing does not depend on evolutionary status of the outflows although 
further studies are required to confirm this. 

• Estimates of the atomic and molecular mass flux rates have been derived from the 
luminosity of the outflows and the kinematical information given by sub-mm interfer- 
ometric data. Values of the order of ~ 10~^ and ~10~^ Mq yr~^ are inferred from the 
[SI]25//m and H2 line luminosities, respectively. A comparison with the momentum 
flux of the large scale CO molecular outflow suggests that the detected atomic jet has 
the power to drive the large scale outflow and thus may represent the primary jet 
ejected by the source. 
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Actions - Human Resource and Mobility within the JETSET (Jet Simulations, Experiments 
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Table 1. L1448 observed lines:SPITZER/IRS 
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'^the [SI]25.2/xm line is detected in two different orders of the LH module. 
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Table 2. L1448 observed lmes:UKIRT/UIST 
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s 


(0) 


< 0.4 


0.35±0.05 


2 


10±0.06 


2.2477 


H2 


2-1 


s 


(1) 


< 0.4 


< 0.05 





97±0.04 


2.3556 


H2 


2-1 


s 


(0) 


< 0.3 


< 0.04 





32±0.04 


2.4066 


H2 


1-0 


Q 


(1) 


8.07±0.39 


1.31±0.04 


9 


73±0.04 


2.4134 


H2 


1-0 


Q 


(2) 


3.28±0.36 


0.45±0.05 


4 


10±0.03 


2.4237 


H2 


1-0 


Q 


(3) 


8.01±0.37 


1.56±0.05 


9 


86±0.04 


2.4375 


H2 


1-0 


Q 


(4) 


1.49±0.28 


0.42±0.05 


2 


59±0.06 


2.4548 


H2 


1-0 


Q 


(5) 


5.89±0.30 


1.09±0.05 


5 


34±0.04 


2.4756 


H2 


1-0 


Q 


(6) 


< 0.4 


< 0.05 


1 


75±0.05 


2.5001 


H2 


1-0 


Q 


(7) 


< 0.5 


< 0.06 


3 


56±0.06 



Table 3. Temperatures and Column Densities 





H2 




Fell / Sill 


Position 


T(K) 


N(H2) (cm-2) 


T(K) ne(cm-3) 


CS 


900±170 


3.18 10^^ 




CN 


755±35 


7.81 10^^ 


<1700 200-1000 


OFl 


600±30 


1.92 10^8 




0F2 


670±50 


2.34 lO^s 


<2500 300-500 



-23- 



Table 4. Gas-phase abundance of Fe and Si 



Position CN OFl 0F2 

[Fegas\/[Fe(:X 3-5 IQ-^ 4-7 IQ-^ 3-5 IQ-^ 

[Fegas]/[FeQ]^ 1-2 IQ-^ 1.5-3 IQ-^ 1-2 10"^ 

[Sigas]/[SieT 3-510-2 ■■• 3-5 10"^ 

[SigasMiSiet 1-210-1 •■• 1-210-1 

" assuming collisions with electrons and Ug =400 cm"^ 
^ assuming collisions with n^, njot^lO^ cm-^ and tlh — 0.1 n^t 



Tabic 5. Physical parameters 



Position 


M(H2) 


M([SI]) 


M([FeII]) 






(lO-^Moyr- 




CN 


0.5 


9-20 


2-20 


0F2 


1.4 


3-8 


0.9-9 




Fig. 1. — Positioning of IRS slits for the SL, SH, & LH modules, overlayed on 4.5 /im 
IRAC image (from the c2d survey). For the SL and LH modules, both nodding positions 
are displayed, whereas for the SH module only the outline of the observed area is presented 
for clarity. The i ntersections of the slit s define four areas of extraction along the jet axis. 
(Upper left) SiO (iGuilloteau et al.lll992l ) contours overlayed on the 4.5 /im image show that 
the driving source is almost coincident with the northern emission area; crosses indicate the 
centers of the CS, CN and OFl positions from bottom to top. 
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Fig. 2. — Extracted spectra (lower panel: 5-18 fim, upper panel: 18 - 38 fim) for each 
of the four positions defined in Figure [H after defringing and removing zodiacal hght. The 
observed hues are labeled on the top of each set. Open squa res and triangles are I RAQ 
bands 3,4 and MIP S band 1 continuum flux measuremens from |j0rgensen et al.l (120061 ) and 
Tobin et all J2007h . 



- 26 - 




1 I I I I I I I I I I I I I I I I I I I I I I I I I I I r 

• [FeII]26/zm/H2 S(l) 
8 h o [Sl]26/zm/H2 S(l) 

,CN 



.2 6 h 
-i-> 

CO 

u 

c 

J 4 h 



2 - 



I 



OFl 





-20 -10 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 



0F2 

CD 

1 



10 20 

Distance (arcsec) 



30 



Fig. 4. — Ratio of the atomic ([FeII]26/im and [SI]25yum) over the molecular (H2 S(l)) 
line emission. Atomic emission dominates close to the driving source (CN) while the 
atomic/molecular contribution decreases in the OFl and 0F2 positions (see Fig. [1]). 
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Fig. 5. — Excitation diagram of the H2 lines for each position defined in Fig. [H Spitzer (open 
circles) and UKIRT (open squares and arrows for the upper limits) data points are optimally 
least square fitted (solid line). Temperatures and total column densities are derived from 
the slope of the fitted line, and its intersection with the ln(Ny^j/g) axis. Visual extinction 
values used for dereddening and derived temperatures are displayed in the upper left corner 
of each panel. 
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Fig. 6. — Column Density Ratio (CDR) diagram of the 0F2 position, constructed normal- 
ising the column density of each transition to the value given by a gas at lOOOK and to 
the 1-0 S(l) column density. Triangles represent Spizer 0-0 lines, squares represent lines 
from the first vibrational levels, including the 1-1 upper limits derived from the Spitzer IRS 
spectra, and diamonds correspond to 2-1 S lines. The superimposed model predictions are 
for the ground (sohd), first (dashed) and second (dot-dashed line) vibrational level. The fit 
corresponds to a bow shock moving at 100 km s~^ into a medium of total density 10^ cm~^. 
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Fig. 7. — Diagnostic diagram of the [SiII]35yum/[FeII]26/im ratio versus electron density, for 
temperatures of lOOOK and 2000K (solid and dashed lines). Hatched areas represent the 
observed ratio for the CN and 0F2 regions that correspond to electron densities of from 200 
to 1000 cm~^ and from 300 to 500 cm~^ respectively. 
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Fig. 8. — Diagnostic diagram of [FeII]18/im/[FeII]26/im ratio versus temperature, for electron 
densities of 2 10^, 5 10^ and lO^cm"^ (solid and dashed lines respectively). The non-detection 
of the[FeII]18/im line allows us only to constrain the upper temperature limits for the CN 
and 0F2 positions 



